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Abstract 

The condition for pressure isotropy, for spherically symmetric gravitational fields with 
charged and uncharged matter, is reduced to a recurrence equation with variable, rational 
coefficients. This difference equation is solved in general using mathematical induction 
leading to an exact solution to the Einstein field equations which extends the isotropic 
model of John and Maharaj. The metric functions, energy density and pressure are well 
behaved which suggests that this model could be used to describe a relativistic sphere. 
The model admits a barotropic equation of state which approximates a polytrope close 
to the stellar centre. 
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1 Introduction 



Solutions of the Einstein field equations for spherically symmetric gravitational fields in static 
manifolds are necessary in the description of compact objects in relativistic astrophysics. The 
models generated are utilised to describe relativistic compact objects where the gravitational 
field is strong as is the case in neutron stars. It is for this reason that considerable energy and 
time is devoted to the study of the mathematical properties and features of the underlying non- 
linear differential equations. The detailed lists of Stephani et al PQ and Delgaty and Lake [2] for 
static, spherically symmetric models provide a comprehensive collection of interior spacetimes 
that match to the Schwarzschild exterior spacetime. It is important to note that only a few of 
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these solutions correspond to nonsingular metric functions with a physically acceptable energy 
momentum tensor. Some of the exact solutions to the field equations, which satisfy all the 
physical requirements for a relativistic star, are contained in the models of the Durgapal and 
Bannerji [3], Durgapal and Fuloria [I], Finch and Skea jS], Ivanov [6 , Maharaj and Leach [7j 
and Sharma and Mukherjee [S], amongst others. 

In this paper our objective is to find new exact solutions to the Einstein field equations 
which may be used to describe the interior spacetime of a relativistic sphere. The approach 
essentially reduces to the analysis of difference equations which we demonstrate leads to explicit 
solutions. We first express the Einstein equations, for neutral matter, as a new set of differential 
equations utilising a transformation due to Durgapal and Bannerji |3] in §2. We choose a 
general polynomial form for one of the gravitational potentials, which we believe has not been 
studied before. This enables us to simplify the condition of pressure isotropy in §3 to a second 
order linear equation in the remaining gravitational potential. We assume a series form for 
this function which yields a difference equation which we manage to solve using mathematical 
induction. It is then possible to exhibit a new exact solution to the Einstein field equations 
which can be written explicitly as shown in §4. Our results contain the model of John and 
Maharaj [9] as a special case. The curvature and matter variables appear to be well behaved in 
the interior spacetime. We also demonstrate the existence of an explicit barotropic equation of 
state relating the pressure to the energy density. For small values of the radial coordinate, close 
to the stellar core, we demonstrate that the equation of state approximates a polytrope. Then 
in §5 we consider the Einstein-Maxwell equations for charged matter. The solutions found in 
the presence of an electromagnetic field reduce to the model for neutral matter given earlier in 
§4. We believe that the method illustrated in this paper is a useful device in the production of 
models for compact objects. 



We assume that the spacetime manifold is static and spherically symmetric. This requirement 
is consistent with models utilised to study physical processes in relativistic astrophysical objects 
such as dense stars. The generic line element for static, spherically symmetric spacetimes is 
given by 



2 Field equations 



ds 2 = - e Mr) dt 2 + e 2X(r) dr 2 + ^2 + ^2 g #2) 
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in Schwarzschild coordinates. For neutral perfect fluids the Einstein field equations can be 
expressed as follows 
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for the spherically symmetric line element The energy density p and the pressure p are 
measured relative to the comoving fluid 4-velocity u a = e~ v b\ and primes denote differentiation 
with respect to the radial coordinate r. In the field equations (0) we are utilising units where 
the coupling constant = 1 and the speed of light c = 1. 

A different but equivalent form of the field equations is obtained if we introduce a new 
independent variable x, and new metric functions y and Z, as follows 

A 2 y 2 (x) = e 2u{r \ Z(x) = e~ 2X{r \ x = Cr 2 (3) 

In the transformation (JHJ), due to Durgapal and Bannerji the quantities A and C are 
arbitrary constants. Under the transformation (|3*|). the system (J2J) has the form 

1^-2Z = £ ,4a) 

4Z» + ^i = £ (4b) 
y x C 

AZx 2 y + 2Zx 2 y + (Zx- Z + l)y = (4c) 

where the dots denotes differentiation with respect to the variable x. It is clear that (plj is 
a system of three equations in the four unknowns p,p,y and Z. The advantage of using (Jlj). 
rather than the original system (J2J), is that on a suitable specification of Z(x), (|4*cj) becomes 
second order and linear in y. 



3 Difference equations 

The approach that we follow is to specify the gravitational potential Z[x) and attempt to solve 
()4c|) for the potential y. From inspection it is clear that the simplest solutions to the system 
P|) correspond to polynomials forms for Z(x). Consequently in an attempt to obtain a new 
solution to the system (JH) we make the particular choice 

Z = 1 + ax k (5) 
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where a is a constant. As far as we are aware all exact solutions found previously correspond 
to forms of the gravitational potential Z(x) which are linear (k = 1) or quadratic (k = 2) 
in the independent variable x. For a comprehensive list of spherically symmetric solutions 
corresponding to different choices of Z(x), relevant to relativistic astrophysics, see Delgaty 
and Lake [2] and Stephani et al pQ. Recently John and Maharaj |U] generated a new solution 
corresponding to the cubic form with k = 3. Higher values for k have not been considered before 
because the resulting differential equation in the dependent variable y is difficult to solve. With 
the specified function Z given above the condition of pressure isotropy (|4"c|) becomes 

2 (l + ax k ) y + akx k ~ l y + aax k ~ 2 y = (6) 

where we have set a = for convenience. The linear second order differential equation (jHJ) 
is difficult to analyse when a ^ 0. Standard handbooks of differential equations and computer 
software packages such as Mathematica have not been helpful in generating exact solutions to 
(0). We attempt to find a series solution to (jUJ) using the method of Frobenius. As the point 
x = is a regular point (for k > 2) of ©, there exist two linearly independent solutions of the 
form of a power series with centre x = 0. Therefore it is possible to write 

oo 

y{x) = c nX n (7) 

n=0 

where the c n are the coefficients of the series. For an acceptable solution we need to determine 
the coefficients c n explicitly. 

Substituting ((7|) into ^ yields 

fe-i 

2 Y c nn{n - l)x n - 2 + [2c k k{k - 1) + aac }x k - 2 

n=2 

+ [2c k+1 (k + l)k + ac x {k + a)]x k ~ l 

oo 

+ ^2[2c n+k (n + k)(n + k-l) + ac n (2n(n-l) + kn + a)}x n+k ~ 2 = 0. 

n=2 

Note that in this equation each term contains ascending powers of x. For this result to hold 
true for all x we require that the coefficients satisfy a set of consistency conditions 

c n = 0, n = 2, 3, .., k — 1 (8a) 

2c k k(k - 1) + aac = (8b) 

2c k+1 (k + l)k + aci(k + a) = (8c) 

2c n+k (n + k)(n + k - 1) + ac n [2n(n - 1) + kn + a] = 0, n > 2. (8d) 
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The recurrence relation (|8d|) is of order k and consists of variable, rational coefficents. Equation 
(|8dp does not fall in the known classes of difference equations and needs to be analysed from 
first principles. Note that ()8aj) and (|8djl imply 

Cfc+2 = Ck+3 = • ■ ■ = C fc+ (fc_i) = (9) 



Hence the remaining nonzero terms are applicable from n > k in (j8d|) . Consequently the system 
((HI) is reduced to the following set 

a a 

Ck = -2k(k^) C0 ( } 

a (k + a) , , , 

c « = "2liTT)I c ' < 10b > 



2n(n — 1) + kn + a 



c n , n > k. (10c) 



_ (n + fc)(n + A; - 1) 

Note that from (J%a|) and © we have two sets of consecutive coefficients which vanish. This 
pattern of zero coefficients repeats itself because of the recurrence relation (jlUcJl . The nonvan- 
ishing coefficients can be placed into two groups. These can be written either in terms of the 
first leading coefficient cq or in terms of the second leading coefficient c%. We consider these 
sets in turn. 

We first consider the coefficients 

Co, Cfcj C2k, C3fc; ■ • • 

From the system (fTUj) we can generate expressions for c&, C2k, c^, ■ ■ ■ m terms of the first leading 
coefficient Cq. These coefficients generate a pattern and we can write 

/ ay+^ [2(kp)(kp-l) + k(kp) + a} 
Cnk+k = ("2 ) l\ (k P + k)(k P + k-l) C °' H ~ ° 

where we have utilised the conventional symbol Yl f° r multiplication. It is also possible to 
establish the result (fTTj) rigorously applying mathematical induction. For n = the result (jllj) 
is obvious since 

a\0+i 2.0(0- l) + k.0 + a 
\) 



n 1 2) (o + jfe)(0 + fc-l) C ° 



Now suppose that the result (jlljl holds for n = r which is the inductive step so that we can 
write 

( a V +1 TT t 2 (M(^- 1) + K k P) + «] 
c rk+k - (--j 11 _______ Co . 

Then let n = r + 1 in (|11|) which is the next term. From equation ()10c|) we have that 



a [2k(r + l)(k(r + 1) - 1) + k(k(r + 1)) + a] 

C {r+ l) k+k - ~- (fc( r . + l) + j t )(A;( r . + l) + jfc_l) ° rk+k 
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By the above inductive step this is equivalent to 

a [2k(r + l)(k(r + 1) - 1) + k(k(r + 1)) + a] 
c [r+1)k+k - -- (k( r + 1) + k)(k(r + 1) + k - 1) X 

ay+ l [2(kp)(kp - 1) + k(kp) + a] 

11 ~ ~^TT, — ~ C 



u (kp + k)(kp + k - 1) 



ay+ 2 j± [2(kp)(kp - 1) + k(kp) + a] 
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The above equation shows that the result holds for n = r + 1. Hence by the principle of 
mathematical induction the result (JTTj) is true for all nonnegative integers n. 
We can perform a similar analysis for the coefficients 

Cl, Cfc + i, C2fc+l, C3fc + i, . . . 

and attempt to establish a general pattern. From the system (fTU|) we can obtain the expressions 
for Cfc+i, C2fc+i, C3fc+i, . . . in terms of the second leading coefficient ci. These coefficients also 
generate a similar pattern as above and we can write 

a V+ i " [ 2 (fcp + 1)(M + + 1) + a] 

11 (]fcp + ]fe+l)(Jfcp+Jfc) Cl ' ^ 



and n denotes multiplication. As in the previous case it is possible to verify the result (fT2|) 
rigorously. For n = the result (|12|) is obvious since 

/ a\0+i [2(jfe.O + l)(fc.O) + fc(fc.O + 1) + a) 
c k+ i - [--) (jfc.O + fc + l)(A;.0 + A;) ° X 

Now suppose that the result (|12j) holds for n = r. Then we have 

a V +1 TT [ 2 ( k P +1 )( k P) + H k P + l) + a] 
2 J 1 = 1 (A;p + A; + l)(fcp + ifc) Cl ' 

Now let n = r + 1 in (|12|) which is the next term. From equation (jlUcj) we have 

a 2[(r + l)k + l](r + l)Jfc + fc[(r + l)fc + 1] + a 

C (r+1)fc+fe+1 - -- [k{r + 1) + k + lMr+1) + k] C (r + Dk + 1- 

By the above inductive step this is equivalent to 

a 2[(r + l)k + l](r + l)k + fc[(r + l)Jfe + 1] + a 
c (r+1)k+k+1 - -- [A;(r + 1) + A . + 1][A . (r + 1) + A;] x 
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The above equation shows that the result ()12|) holds for n — r + 1. Hence by the principle of 
mathematical induction the result (|12j) holds for all nonnegative integers n. 

Thus we have solved the general recurrence relation in (JHJ) and The coefficients 

Cfc, c 2 fc, c 3 fc, . . . are generated from (fTTj) . The coefficients Ck+i, c 2 k+i, C3k+i, ■ ■ • are generated 
from (|T2*|) . Hence the difference equation (jlUcjO has been solved and all nonzero coefficients are 
expressible in terms of the leading coefficients Cq and c±. From the equations (J7J), (|TT|) and ()12j) 
we can write the general solution to © as 
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where c and ci are arbitrary constants and a = k -^-. Clearly the solution (fTHJl is of the form 
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(14) 



(15a) 



(15b) 



are linearly independent functions. Therefore we have found the general solution to the differ- 
ential equation (JBJ) for the particular gravitational potential Z given in 

From (|14j) and (|15|) we can generate a number of new particular exact solutions for specified 
values of k and a. For certain values of k the solution may reduce to models that have already 
been found. We let k = 3 so that a — 1. Then the gravitational potential Z becomes 



and the gravitational potential y is 

oo 



Co 



Ci 



;r 



n=0 

oo 

+E 
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ax 
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(16) 

Therefore our general solution contains as a special case pfijl which is the model of John and 
Maharaj jHj- 
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4 Einstein models 



From the analytic representation (|14|). (fTo")) and the Einstein field equations (jlj) we generate 
the exact solution 



e 2X = i 1 , (17a) 

1 + ax k v ' 

e 2v = A 2 y 2 (17b) 

-a(l + 2k)x k ~ 1 (17c) 



£ = A(l + ax^ + ax^ 1 (17d) 

where A and C are constants, y is given by (fl"4"|) . a is constant and A; is a natural number 
greater than one. This solution has a simple form and is expressed completely in terms of 
elementary functions and the series (fTKj) . The metric functions and the matter variables are 
written explicitly in terms of the independent variable x. 

The gravitational potentials v and A satisfy the conditions of physical reasonability for a 
relativistic sphere: they are finite at the centre x = 0, continuous in the interior of the spacetime 
and match smoothly to the Schwarzschild exterior spacetime at the boundary of the sphere. To 
obtain positive energy density we require a < 0. The matter variables p and p are bounded and 
nonsingular at the origin. Even though our solution is presented in terms of an infinite series, 
we point out that by using computer software packages it is easy to generate graphical plots of 
v, A, p and p. Thus a physical analysis of our model is very feasible and this is an avenue for 
future research. The expressions given above for the gravitational potentials and the matter 
variables have the advantage of simplifying the analysis of the physical features of the solution, 
and will assist in the description of a relativistic compact bodies such as neutron stars. 

From ()17|) we can observe that 

x = I — I , a < 0. 



-aC(l + 2k) 

Therefore from equation (j!7d|) . we can express p in terms of p only. This is a special feature 
of our model and most of the solutions presented in the literature do not satisfy this property. 
As this property does not depend on a particular value of k, it is valid for allowable values of 
k. Hence the solution satisfies the barotropic equation of state 

p = p(p). 

This property is usually made as a requirement for a physically relevant model. Also observe 
that for small values of x close to the stellar centre we have y ~ Cq + C\X. Then from (|T7|) we 
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have the approximation 

L « *2± (18) 

C + Ci 



-oC7(l+2fc) 

Therefore for small values of x close to stellar centre ()18|) implies that we have the approximate 
equation of state 

V oc 

which is of the form of a polytrope. We make the observation that equations of state, with 
pressure as a negative power of energy density, arise in relativistic cosmological models with 
a Chaplygin gas; these models are now widely used as a driver of acceleration to generate the 
present universe [THl HI] . 



5 Einstein- Maxwell models 

It is possible to extend the solutions presented in this paper to include the electromagnetic field 
as the procedure to obtain the new general solution is very similar to §3. We present only an 
outline and do not give all the details of the calculation. A generalisation of the Einstein field 
equations is given by 

1 - Z ■ p E 2 , . 

—~ 2Z = h + 2d (19a) 

iZ i + ^i = i-i; (18b) 

y x C 2C 

4Zx 2 y + 2Zx 2 y+lZx-Z + l-—jy = (19c) 

a 2 AZ f ■ \ 2 

- = -(xE + B) (19d) 

where E is the electric field intensity and a is the charge density. When the electric field E — 
then the Einstein- Maxwell equations (fTU|) reduce to the Einstein equations for neutral 
matter. The system of equations ()19j1 governs the behaviour of the gravitational field for a 
charged perfect fluid. 

On substituting (J3J) in equation (|19c|) we obtain 

2 (1 + ax k ) y + akx k - l y + (^^x k ~ 2 - y = 0. (20) 

If we specify the electric field intensity to be 

E 2 = 2aC(5x k ~ l 
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where (3 is a constant, then (}2T)j) becomes 



2(1 + ax k ) y + akx k 1 ij + a(a - (3)x 



k-2„ 



(21) 



where a = ^p. This choice for E is physically reasonable as it is finite at the centre and 
continuous in the interior spacetime. It is now possible to use the results of §3 to solve ()21j) 
directly. On comparing the differential equations (JHJ) and (}2T|) . we can write the solution to the 
differential equation (|2"T|) as 

y{x) = c yi(x) + ciy 2 (x) (22) 



where 
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(23b) 



From the analytic representation 
generate the new exact solution 
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and the Einstein-Maxwell field equations f)19|) we 



(24a) 

(24b) 
(24c) 

(24d) 



1 + ax k 
AY 

-a(l + 2k + (3)x k ~ 1 

4(1 + ax k )- + a(l + (3)x k ~ l 

y 



— = 2a(3x 



k-l 



(24e) 



where A, C, a, and (3 are constants, y is given by (|2~2~|) and A; is a natural number greater than 
one. For positive energy density we require a < 0. We note that we can express p in terms of p 
only; hence this model also satisfies the barotropic equation of state p = p(p) as in the case of 
neutral fluids. When (3 = 0, reduces to (|T7|) which is the uncharged solution found earlier. 
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